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ABSTRACT 

In this paper we study the cosmological consequences of the minimal su- 
persymmetric extension of the standard model in the case that the neutralino is 
heavier than the W . We calculate the cross section for annihilation of heavy neu- 
tralinos into final states containing gauge and Higgs bosons (xx WW, ZZ, 
HH, HW, HZ) and combine these results with those previously obtained for 
annihilation into fermions to find the relic cosmological abimdance for the most 
general neutredino. The new channels are particularly important for Higgsino-like 
and mixed-state neutralinos, but are sub-dominant (to the fermion-antifermion 
ann ihilation channels) in the case that the neutralino is mostly a gaugino. The 
effect of the top quark mass is also considered. Using these cross sections and 
the cosmological constraint < 1, we map the entire range of cosmologically 

acceptable supersymmetric parameter space and discover a very general bound 
on the neutralino mass. For a top quark m 2 iss of less than 180 GeV, neutralinos 
heavier than 3200 GeV are cosmologically inconsistent, and if the top quark mass 
is less them 120 GeV, the bound is lowered to 2600 GeV. Neutralino states that 
axe mostly gaugino are constrained to be lighter than 550 GeV. We find that a 
“heavy” neutralino (m^ > mw) that contributes ~ 1 arises for a very wide 
range of model parameters and makes, therefore, a very natural and attractive 
dark matter candidate. 
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I. INTRODUCTION 


Low-energy-supersymmetric theories provide an elegant solution to the hi- 
erarchy problem and have been extensively studied in recent years. In most 
models, the lightest supersymmetric particle (LSP) is stable and makes an at- 
tractive candidate for the dark matter known to exist throughout the Universe 
and especially that in galactic halos?’^ In the minimal supersymmetric model, 
the most likely LSP is the neutralino— a linear combination of photino, Higgsino, 
and Z-ino fields. In addition, it has been shown that if it exists and is stable, 
it very naturally has a relic abundance near closure density, and thus proba- 
bly comprises a significant component of our galactic halo independent of the 
parameters of the supersymmetric modelf’^ However, in previous work it has 
been implicitly assumed that the mass of the neutralino, rrijf , is less than the W 
mass. (But see the recent work of Olive and Srednidci.® ) Important new anni- 
hilation channels open up when > mw and because the relic abundance is 
determined by the total annihilation cross section, it is crucial to include these 
new channels. In particular, the results of Refs. 4 and 5 were found considering 
only the reaction XX fL where / is a quark or lepton, while for > mw 
final states such as W^W~ , ZZ, and states involving Higgs bosons must also be 
included. For Dirac neutrino dark matter, it was shown by Enqvist, Kainulainen, 
and Meialampi^ that the gauge-boson final states dominate as the neutrino mass 
increases, and we show that the same is true for the neutralino in some regions 
of supersymmetric parameter space. 

Since a supersymmetric solution to the hierarchy problem requires a relatively 
light LSP, an “extremely” massive neutralino is not very attractive; however, 
masses in the 100 GeV range are entirely reasonable, especially as unsuccessful 
accelerator searches push the possible masses of the supersymmetric particles 
upward. In this paper, we consider in detail the possibility of neutralino dark 
matter heavier than the W, t.e., > mw- Regions of parameter space that 

are ruled out due to an excessive relic density of neutralinos are identified and 
the cosmologically viable range of neutralino masses is found to be in a range 
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that is also suitable for solving the naturalness problem. We also show that the 
cosmologiceilly acceptable supersymmetric models naturally predict relic abun- 
dances in the range appropriate for supplying the bulk of the dark matter. We 
find both of these results encouraging. 

The plan of this paper is as follows; In Section II, we describe the low-energy- 
supersymmetric model upon which our calculations axe based and describe the 
calculation of the relevant annihilation cross sections. The complete calcula- 
tion is done for the W~^W~, ZZ, and neutral Higgs final states, and reasonable 
estimates are made for the charged Higgs and mixed Higgs-gauge-boson final 
states. We then compare the cross sections for these channels to those previously 
computed for the fermion-antifermion final states and map out the regions of 
supersymmetric parameter space where the new channels are dominant and the 
regions where they can safely be ignored. In Section III, we identify those regions 
of M-fi space that are ruled out by cosmological considerations and obtain a max- 
imum mass for a cosmologically viable neutralino. In Section IV, the rem8iining 
cosmologically acceptable regions of parameter space are mapped out and it is 
shown that for most models, an interesting relic abimdance (0.01 < < 1) 

is possible, if not probable. Section V is a summary of the paper and the four 
Appendixes contain the detailed expressions for the cross sections as well as their 
forms in several useful limiting cases. 

II. MODEL AND CROSS SECTIONS 

We start with the minimal supersymmetric extension of the standard model 
as described by Haber and Kane/ and follow their conventions and notation 
throughout. In these models there are four neutralinos, linear combinations of 
the piirtners of the photon, Z, and two neutral Higgs bosons; we refer to the 
lightest (the neutralino) as the neutralino and denote it as x, 

x-=Zn,S + z„iW^ + z„3Sj + z,A, ( 1 ) 

where (Z),j is the real orthogonal matrix that diagonalizes the neutrahno mass 
matrix. In addition, there are two chaxginos, linear combinations of the partners 
of the charged Higgs and gauge bosons. The neutralino and chargino masses and 
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mixings are determined by three^^ arbitrary parameters: two mass parameters, /x 
and M, and the ratio of Higgs vacuum expectation values, tan/0 = V 2 /v\. In ad- 
dition, there are scalar quarks and scalar leptons whose masses are arbitrary, but 
to simplify matters, a single common squark and slepton mass rrip is assumed. 
Finally, there are the two requisite Higgs doublets which give rise to two neutral 
scalars and if®’ ^ neutral pseudoscalar and two charged Higgs bosons, 
H'^ and H~. The masses and couplings of these particles are determined by /z, 
M, and tan as well as the mass of the lightest scalar mjjo (see Refs. 11 and 12). 
In these models the masses of the neutral Higgs particles satisfy: mjfo > mz] 
mjjo < mz cos 2^; and mgo < m^o < m^o. The charged Higgs particles are 
both heavier than the W. Feynman rules as well as mass matrices and other 
relevant formulas for the Higgs particles can be foimd in Refs. 11 and 12. 

Although the analysis of the heavy neutralino can become very involved, 
many results can be understood in terms of the neutralino properties displayed 
in Fig. 1. The broken cvuwes in Fig. 1 are neutralino mass contours in the M- 
/i plane, and the solid curves are the gaugino fraction (defined as + ZI 2 ) 
of the neutralino. In Fig. 1, tan^S = 2 is assumed; however, these curves are 
relatively insensitive to the value of tan;d. The mass contours are hyperbolas 
that asymptote to = |/x| for large M and to = M' th M/2 for large ||x|. 
For large values of i/x| and M, models where the neutralino is half Higgsino and 
half gaugino fall along the line /x = \M tan^ $w « M/2. In the regions where 
the gaugino fraction is greater than 0.99, \ is almost a pure R-ino state {Zni ^ 1 
and Zni ^ 0 for t ^ 1). (When M' = |Mtan^ Bw, the lightest neutralino cannot 
both be a pure photino and heavier than the W.) When the gaugino fraction is 
less than 0.01, \ is very nearly a pure Higgsino state (Z„j « Z„2 ^ 0). 

In order to csilculate the relic abimdance we need the total cross section 
for annihilation of neutralinos into all lighter particles. The result for // final 
states, where / is a quark or lepton, has been calculated previously* and appears 
(with corrections) in Appendix D. Other possible final states include the gauge 
bosons, W'^W~ and ZZ, the six possible combinations of two neutral Higgs 
bosons H\H 2 , etc.), the charged Higgs, and the five 
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Higgs-gauge-boson final states, W'^H , W ZHi, ZH 2 , and ZH^- 


Consider the gauge-boson final states first. The Feynman diagrams for \X 
W'^W~ are shown in Fig. 2, and the diagrams for XX are shown in Fig. 3. 

There are a total of seven diagrams contributing to the WW final state and ten 
to the ZZ final state. The complete cross sections are given in Appendixes A 
and B. The calculation of these cross sections is a tedious process with ample 
opportunity for error; moreover, the final expressions for the cross sections are 
long. Fortimately, unitarity provides a very nice check. The contributions of 
the individual Feynman diagrams in Fig. 2 to the cross section for annihilation 
into W'*^W~ contain terms proportional to s, the center of mass energy squared. 
However, each partial wave of the complete matrix element must be bounded 
by a constant which implies that as s — > 00 , the total cross section must be 
proportionjil to For this to occur, a highly nontrivial cancellation must 

take place between the squares of the seven individual diagrams and the twenty- 
one interference terms. The formulas for the total cross section presented in 
Appendix A (as well as those for the ZZ final state presented in Appendix B) 
do indeed exhibit this behavior. In addition, in the low-ener^ limit, relative 
velocity v — » 0, the total cross sections reduce to simple squares, Eqs. (A 10) and 
(B7). The correct high-energy behavior of the total cross sections as well as the 
fact that the cross sections simplify in the u -♦ 0 limit give us some confidence 
that our complicated expressions are correct. 

The Feynman diagrams for xX are shown in Figs. 4 and 5, and 

the resulting cross sections are given in Appendix C. None of the squares of the 
individual diagrams for these channels diverge for large s, so unitarity cannot 
be used to check our results. Moreover, there is no particular use in obtaining 
expressions that are exact in s. Therefore, we have obtained expressions that 
are correct only to order which simplifies matters greatly. The Higgs-boson 
channels tend to be sub-dominant to the gauge-boson and fermion-antifermion 
rhaTtupls for most of the heavy neutralino parameter space; however, some of 
the HfHj final states are quite important when rrij^ < mw and the effect of 
Higgs-boson final states on light neutralinos is considered in another paper As 
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a check, a symbolic manipulation program (MACSYMA) was also used to derive 
these cross sections and numerical agreement between the MACSYMA results 
and the formulas of Appendix C was foimd. 

We have completed preliminary calculations for the mixed Higgs-gauge-boson 
final states, but do not include those results here because these chamnels seem to 
be sub-dominant. Since the matrix elements for these channels contain couplings 
that appear in the pure-Higgs or gauge-bc«on final states, the cross sections for 
the mixed channels should lie in between those for the Higgs-boson final states 
and those for gauge-boson final states. Since the cross sections for the Higgs- 
boson final states and those for the gauge-boson final states are rarely of the 
same magnitude for heavy neutralinos, the mixed states should be negligible. 
Numerically, our preliminary results support this hypothesis. 

To .see the importance of the bosonic final states, we compare the newly 
calculated cross sections with the previously calculated annihilation cross section 
(// final states only).^ Unfort imately, even with the simplifications already made, 
a six-dimensional parameter space (^, M, tan/3, m^o, and the top quark 
mass, mt) must be explored. To simplify matters we consider only three values of 
rrij: my = m^, my = 2m^, and my = oo. Recall that m^j provide a lower limit 
to the squark mass since x is assumed to be the LSP. When using my = mj, the 
importance of the new (bosonic) final states will be underestimated because m y 
appears only in the denominator of the cross section into //. The opposite limit, 
my = oo, supresses the fermion channels greatly, so by considering both limits 
the sensitivity of our results on m y/ m^j can be seen. 

To simplify further, we will set m^o = 0. The lightest Higgs mass mffo 
determines a Higgs mix ing angle a, and coupling constants of the Higgs bosons 
depend only trigonometrically on a. When annihilation into all bosonic final 
states is considered, the total cross section is relatively insensitive to the choice 
of mgo provided that the neutralino is heavy enough that threshold effects are 
unimportant. (As discussed in Ref. 9, the value of rngv is important for light 
neutralinos.) For the cases presented below, we have confirmed numerically that 
the resxdts are insensitive to value of mjyo. We have included Higgs exchange in 
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the fermion channels, but for very massive neutralinos this effect is not important. 
As discussed in Ref. 9, Higgs exchange can be important, but only when neeir a 
pole. 

Using = ruj^ and mjyo =0, we show in Fig. 6 the regions of the M-n 
plane where the new final states are important for several values of tan/? and 
mi and for a neutralino-neutralino relative velocity of « 1/4 (relevant for 
annihilation in the early Universe around “freeze-out”). In the regions “hatched” 
with positively sloped lines, the annihilation cross section into fermion final states 
is greater than that into gauge and Higgs bosons by at least a factor of ten. In 
the areas hatched with negatively sloped hnes, the gauge-boson contribution 
to the annihilation cross section is at least ten times as great as that of the 
Higgs and fermion channels, and in the regions hatched with horizontal lines, 
the Higgs-boson contribution dominates that of fermions and gauge bosons by at 
least a factor of ten. In the areas hatched with vertical lines, none of the three 
contributions dominate by more than a factor of ten. 

The analogous plots for a relative velocity w 0 are very similar and so 
are not shown. We remind the reader that there are two distinct energy regimes 
of interest for XX annihilation: « 1/4, typical of the early Universe around 

“freeze-out” when the annihilation cross section determines the relic abundance 
of neutralinos; and ss 10“® « 0, characteristic of the neutralino annihilations 
in the galactic halo or in the body of the Sun or Earth, which may prove to be 
an interesting soiurce for high-energy cosmic rays or neutrinos. Since the rela- 
tive velocity is small in both regimes, an expansion to first order in powers of 

is a good approximation for the total annihilation cross section ov. In the 
early Universe, one might expect the term in av to be smaller by a factor 
of v^/4 fa 1/16, but for Majorana particles such as the neutralino, the v® term 
may be suppressed and the term may actually dominate. For present day 
annihilations with v^/4 < 10“®, the term in the cross section is never impor- 
tant. Even if the term vanished at tree level, one would expect higher-order 
corrections and/or three-body final states to overwhelm the term in this ve- 
locity regime.^^ Given the different possible structures of the cross section in the 
two velocity regimes, it may be surprising that the figure corresponding to Fig. 6 
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with w = 0 (not shown) looks so similar. In regions of large neutralino mass this 
similarity occurs because the cross section for annihilation into gauge bosons and 
top quarks has no s-wave (v° term) suppression. In the regions of small 
the similarity is due to our assumption that m^o =0. The channel in 

particular is eilways open and has no s-wave suppression. If we use larger (and 
more realistic) values for m jjo so that annihilation of light neutralinos into Higgs 
bosons is kinematically forbidden, the behavior of the total annihilation cross 
section in the galactic halo differs from that in the early Universe. We address 
this further in Ref. 9. 

FVom Fig. 6 we see that the gauge-boson final states dominate the annihilation 
cross section when n is small which corresponds to a “Higgsino-like” neutralino 
state. In this case, Zni Zn2 ^ 0 and \Znz \ « \Zn\\ » and the annihilation 

into W^W~ (ZZ) proceeds through the t- and u-channel exchange of the lightest 
chargino (neutralino) with couplings that only depend weakly upon tan /?. For 

= mj, annihilation into fermions proceeds mainly through t- and u-channel 
exchange of the squark, and from the form of the couplings in Appendix D, we see 
that annihilation into the bb {tt) final state increases (decrea^s) with increasing 
tan$. This qualitative behavior is seen in Fig. 6 where the size of the Higgs- 
dominated region decreases with increasing tan^. The effect of the top quark 
shows up dramatically in Fig. 6 where fermion channels dominate along the top 
quark mass threshold when the top is heavy. When M is small, corresponding 
to a “B-ino-like” state, the relevant couplings for XX W~^W~,ZZ are small 
(see Appendixes A and B) and the gauge-boson final states are unimportant. 

For a neutredino that is either gaugino- or Higgsino-like, the Higgs couplings 
are meetly small and those that are non-negligible appear in diagrams where 
heavy virtual particles (heavier neutralinos) are exchanged, so annihilation into 
Higgs final states is relatively small. For a neutralino state that is neither pure 
gaugino nor pure Higgsino (“mixed state”), Higgs channels become more impor- 
tant but om: numerical results suggest that these final states never dominate for 
very heavy (m^f mw) neutralinos. In the “mixed state” regions, several final 
states generally contribute comparable amounts to the total cross section. We 
remind the reader once again that in Fig. 6 the importance of the new channels 
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has been underestimated since we have set my = of course, in general my 
should be greater than m^ — and my could be much greater them rrij^. 

To illustrate the squark-mass dependence of our results, we set my = 2rrij^ 
in Fig. 7. As expected, the regions of domination by fermion final states shrink 
and the regions of Higgs domination and those where no single channel domi- 
nates grow. In Fig. 8, we minimize the effect of the fermion channels by setting 
my = oo, the opposite extreme of my = mjf. Even though fermions can still 
be produced via a Z or Higgs-boson exchange, we see that fermion final states 
hardly ever dominate in this limit. It is also seen that the Higgs final states axe 
more important for B-ino annihilation, while gauge-boson channels dominate for 
Higgsino states. 

As just mentioned, these cross sections can also be used for neutralino anni- 
hilation in the galactic halo or in the body of the Sun or Earth. There eire several 
interesting possibilities for detecting the products from such annihilations. In 
these cases, the limit v — » 0 is appropriate and the formulas simplify consider- 
ably. The cross sections in this limit are displayed in the Appendixes, Eqs. (AlO), 
(B7), (Cll), and (D6). 

III. COSMOLOGICAL CONSTRAINTS 

Using the cross sections discussed in Section II and given in the Appendixes, 
the relic abundance of neutralinos can be found by integrating the Boltzmann 
equation, 

^ + 3Hrij^ = - {(Tv) [n| - , (2) 

where is the actual number density of neutralinos, is the equilibrium num- 
ber density, H is the expansion rate of the Universe, and {erv) is the total annihi- 
lation cross section, thermally averaged and averaged over initial neutralino spins. 
Qualitatively, the number density of neutralinos tracks its equilibrium value imtil 
the annihilation rate T = drops below the expansion rate H (“freeze- 

out”), after which a relic abundance of neutralinos “freezes-in!’ The freeze-out 
temperature Tj depends logarithmically upon (<rtj), but generally occurs for a 
value Xf = mj/Ty « 0(20-30). 



While it is straightforward to integrate Eq. (2) numerically, because of the 
large number of times we must do so, it is far more convenient to use an analytic 
approximation (good to about 5-10%) for the relic abundance instead. In this 
regard we have generally followed the treatment detailed in Ref. 15. 

After freeze-in, the neutralino abundance per comoving volume no longer 
drops rapidly and it is therefore convenient to rewrite the Boltzmann equation 
using the variable Y = n^/s, where s = 2Tr^g*T^/45 is the entropy density and 
counts the total number of effectively relativistic degrees of freedom. Then 
the fraction of criticrJ density contributed by relic neutralinos today is given by 

O ^2 _ Px _ 

" " Pcrit/^2 - (3) 

« 2.82 X 10® Foo(mj^/GeV), 


where the present critical density is Pcrit — 1.05 x 10^ eV cm the present 

Hubble constant is 100 h km sec“^ Mpc“^ with 0.4 < h < 1, and the present 
entropy density is so = 2970 cm 

Expanding the total annihilation cross section in powers of the neutralino 
relative velocity, 

(TV = a + bv^ H , (4) 


it can be shown that an accurate approximation to the neutralino relic abundance 
is given by 

= 0.264<?y ^ 

The freeze-out epoch Xf is determined by 


Xf 


= ln 


0.0764 mpi{a -t- 6&/r/)c(2 -|- c)m^ 


( 6 ) 


which can be solved iteratively to the required precision. Here mp\ = 1.22 x 
10^® GeV is the Planck mass and c is a constant of order unity whose value 
is determined by matching the approximate analytic solutions for x < Xf and 
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X > Xf. In practice, once chooses c to obtain the most accurate approximation 
to Yoo’i we have used the value c = 1/2 which results in a typical accuracy of 
about 5-10% — more than sufficient for our purposes here. 


Figs. 9, 10, and 11 illustrate the dependence of on M and /i for tan /3 = 2 
and rrii = 60 GeV, for several slices through the M-fi plane. The behavior 
of depends on the mass and gaugino fraction of the neutralino and the 

relative importance of the cross sections into the various final states; therefore 
in the following, the reader is encouraged to frequently refer to Figs. 1, 6, and 
8. In Fig. 9, is shown as a function of M for several values of n under 

the two extreme assumptions, rrij = (solid lines) and — oo (broken 
lines). Assuming mj = m^, the relic abundance grows with increasing M (and 
therefore increasing rrij^ (Fig. 1)) while the neutralino is mostly B-ino (Fig. 6) 
and the fermion channels dominate (Fig. 6). This is because ~ ~ 

~ rrij^ ^ due to our assumption that = m^. The neutralino mass 
density reaches a maximum near the value of M where the neutralino is 

an equal mixture of gaugino and Higgsino and drops quickly as the gauge boson 
channels begin to dominate and increase the cross section (Fig. 6). For larger M, 
the neutralino is mostly Higgsino and its mass does not change as M increases 
further (Fig. 1), so levels off. 


The broken curve in Fig. 9 is similar to the solid curve except that we have 

set = oo, greatly suppressing the fermion channels (Fig. 8). For small M, 

the neutralino is J5-ino-like (Fig. 1) and annihilates primarily into Higgs bosons 

(Fig. 8). Since stays nearly constant, we can conclude that the cross section 

in this region is nearly independent of mj. In fact, examination of the cross 

♦ 

section given in Appendix C shows that annihilation is completely dominated 
by the exchange of the heaviest neutralino which is very nearly a pme Higgsino. 
(Exchange of the lighter neutralinos is suppressed due to the lack of a B-ino-B- 
ino-Higgs coupling.) Since the mass of the heaviest neutralino is very nearly equal 
to fi, we find a cross section that is proportional to For large M, the behavior 
is determined again by annihilation into gauge bosons and the suppression of the 
fermion channels due to large mj makes no difference. 
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In Fig. 10, we show for three slices through the M-fi plane with M 

held constant as // increases, again for = mj (solid lines) and my = oo 
(broken lines). For small //, the neutralino is mostly Higgsino (Fig. 1), annihilates 
primarily into gauge bosons (Fig. 6), and increases with ^ (Fig. 1). Again 
~ (av) ^ ~ /rrij^^ ~ Now, however, the switch to a higgsino 

annihilating into fermions results in a jump in (rather than a drop as in 

Fig. 9). Again, for my = m^ (the solid curves) and large values of /x, 
levels off because m^j becomes constant for large //. The broken curves in Fig. 10 
show when my = oo, and we find that at small fx, the behavior is similar 

to the case where my = rrij^. As increases further, the neutralino becomes 
JB-ino-like, and annihilation into Higgs bosons dominates (the fermion channels 
being suppressed by my = oo). The relic abundance increases with fi because 
{av) ~ fi~^ as mentioned earlier. 

Finally, in Fig. 11 we show a slice of the M-fi plane along the line = 
|Mtan^^W' (which for large corresponds to a neutralino that is an equal 
Higgsino-gaugino mixture) for my = m^ (solid line) and my = oo (broken line). 
For small the relic abimdance dips as new annihilation channels become kine- 
matically allowed; however, as seen in Fig. 1, m^j increases monotonically along 
this line when the neutralino becomes heavy, and increases with increasing 

mjj, as it should. As expected from Fig. 6, the relic abimdance for my = oo is 
slightly higher than that for my = since the fermion and gauge-boson con- 
tributions are generally of the same order of magnitude for a neutralino that is 
an equal Higgsino-gaugino mixture. 

The results shown in Figs. 9-11 are qualitatively the same for ^ < 0 and other 
values of tan/9 and m^. 

The age of the Universe at the present epoch (specified by a photon temper- 
ature T = 2.75 K) decreases with increasing values of Qh^. If one insists that 
the age of the Universe today is greater than 10 Gyr and that h > 0.4, the cos- 
mological bound ^ 1 follows. We now use this limit, ^ 1, to 

exclude regions of the M, /x, tan ^ parameter space. As mentioned earlier, rais- 
ing m y lowers the annihilation rate which in turn increases the relic abundance. 
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Therefore we set my = so that regions of peirameter space are excluded 
independently of the unknown squark mass. 

In order to provide realistic limits, we will consider top quark masses in the 
range 60 GeV < mt < 180 GeV. The lower bound follows from unsuccessful 
experimental searches, and the upper bound follows from limits on radiative 
corrrections to s\n^ 6w (Ref. 17). The range of tan^ considered is 1 < tan^S < 
mt/mi,. The reason is that for mt mj, radiative corrections drive tan ^ = V 2 fv\ 
to a value greater than 1. In addition, in many supergravity models, electroweak 
symmetry bre 2 iking only occurs if tan/9 < m(/mj (Ref. 18). 

In Fig. 12, we have hatched out the regions of the M~n plane (for tan/9 = 2) 
that are cosmologically excluded for our limiting cases, mt = 60 GeV and 180 
GeV. (The results for other values of tan /9 are qualitatively similar. The excluded 
regions for /i < 0 are similar to those for // > 0 and so are not shown.) Prom 
Fig. 12 we see that a Higgsino-like neutralino is constrained to be less massive 
than about 3000 GeV; on the other hand, a B-ino-like neutralino must be less 
massive than about 550 GeV. Note that both bounds become more stringent as 
mt decreases. 

To illustrate the dependence of our results on the assumption that my = m^f , 
we show in Fig. 13 the cosmologically excluded regions for my = oo. The region 
of Higgsino-hke parameter space excluded is very similar to the case shown in 
Fig. 12, but many more B-ino-like states are excluded. This is because the 
fermion annihilation channels are (are not) unportant for a B-ino- (Higgsino-) 
like neutralino. 

As the neutralino mass becomes large, the annihilation cross section de- 
creases. Therefore, the constraint < 1 leads to a maximum cosmologi- 

cally acceptable neutralino mass. In order to find the maximum neutralino mass 
consistent with cosmology we must search the parameter space of M, /i, tan^, 
and mt. We did so as follows. For a given mt and tan^, the M-fi plane was 
searched numerically for the largest value of m^ consistent with ^ 1. For 
mt = 60 GeV, the heaviest possible neutralino is mostly Higgsino, and neu- 
tralino annihilation proceeds mainly into gauge bosons. For larger values of m<. 
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the heaviest possible neutralino is nearly half Higgsino and half gaugino, and 
the fermion annihilation channels are of comparable importance. In Fig. 14, the 
maximum cosmologicadly acceptable neutralino mass is shown as a function of 
tan P for several values of m<. As mentioned before, we consider values of tan 
between 1 and FVom Fig. 14 we see that for m* < 180 GeV, a neu- 

tralino mass greater than 3200 GeV is cosmologically unacceptable. If the top 
quark mass should be about 120 GeV, then the neutralino must have a mass less 

than around 2600 GeV. It is interesting that the cosmological window for the 

8 

neutralino mass roughly coincides with that for solving the hierarchy problem. 
These neutralino mass limits would weaken if larger values of m< or more extreme 
values of tan /? were allowed; we believe, however, that we have been generous in 
the regions of parameter space explored. 

Finally, we should mention the status of the “pure” photino and “pure” 
Higgsino states which are frequently considered in the literature as dark matter 
candidates. A pure photino is defined as Zni = cos6\v, Z«2 = sin^W, and 
Z „3 = Zn 4 = 0, while the Higgsino usually co nsi dered has Z„i = Z„2 — 0? 
Zn 3 = sin /3, and Zni = cos 0. When the neutralino is light, an LSP state that 
is a pure photino plus a correction of order (where is the photino 

mass) exists, and is in fact common. Likewise, in the low mass limit there is a 
state that is close to the Higgsino state defined above. Usually however, the above 
relations are assumed to hold without regard to the neutralino mass matrix from 
which they arise. In fact, given the unification assumption, M' = |Mtan^^pv, 
an examination of the mass matrix shows that neither of these states exists as the 
LSP for > mw- If the unification assumption is relaxed, a photino-like state 
of any mass is possible and the relic abundance for these massive photinos can 
be calculated in the same manner as for massive neutralinos. Using my < rriy, 
mt < 180 GeV, and ilh? < 1, one finds a maximum photino mass of around 
600 GeV, a number very similar to our limit for the more realistic J5-ino state. 
For the Higgsino state defined above, the crucial Z-Higgsino-Higgsino coupling 
becomes Z^^^ — Z^^ = cos^ 2^ and such a state annihilates with the cross section 
of a Majorana neutrino times the factor cos^ 20. Tbr large mass however, this 
particular Higgsino solution no longer occurs, and as discussed previously, the 
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actual large mass Higgsinolike solution has Zn 3 ^ i^n4- So we conclude that 
for discussion of supersymmetric dark matter above the W mass, the photino 
and Higgsino states that are usually considered are not so relevant. 

In closing this section, we should again remind the reader that imphcit in 
our calculation of was the assumption that there was no significant entropy 
production after the freeze-out of the rehc abundances. If there were significant 
entropy production — e.g., due to the quark/hadron transition, electroweah sym- 
metry breaking, or relic particle decays — and the entropy per comoving volume 
increased by a factor of 7, then our estimate for would decrease by the 

same factor of 7. On the other hand, our cosmological upper bound to Qj^h? is 
rather generous: Had we assumed that the Universe today must be at least 13 
Gyr old — a very reasonable lower boimd — the cosmological upper limit to 
would fall to about 0.4. Or if we insist that = 1 and that the age of the 
Universe be at least 10 Gyr, then h < 0.7, and ^ 0.5. 

IV. HEAVY NEUTRALINOS AS DARK MATTER CANDIDATES 

It should be noted that for almost any value of M and // not cosmologically 
excluded, a relic neutralino abimdance of significance for the dark matter prob- 
lem, say 0.025 < < 1, can result. These models are interesting because they 

offer an elegant solution to the dark matter problem and predict the existence of 
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particles in our galactic halo which may be detectable. 

In Fig. 15, we show “scatter” plots of vs. mj for models taken from a 
grid of points in the M-/x plane. For fight neutralinos (mj < mw), the relic abvm- 
dances displayed underestimate the expected values because we set rrij = rrij^. 
Since the cross section for annihilation into fermions ~ (m^ + if 

the squark is x times heavier than the neutralino, the relic abundance for fight 
neutralinos is about a factor of -|- 1)^/4 larger than that shown in Fig. 15 (as- 
suming the Higgs channels are negligible). Some indication of this fact can be 
seen in Fig. 16 where we have set = 2m j. As expected, the relic abimdance 
for neutralinos that are mostly B-ino increases roughly as described above. For 
reference, the relic abundances that result when we set = 00 are shown in 
Fig. 17. In Figs. 15, 16, and 17 we arbitrarily set m go =0 since for heavy neu- 
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tralinos, the relic abvmdance is insensitive to this parameter. For light neutralinos 

9 

this is not the case. 

Now consider the relic abundance of very heavy neutralinos. Note that for 
large masses, the relic neutralino abundance seems to cluster along two lines. 
Along the upper line, the neutralino is mostly B-ino. Since annihilation for a 
B-iho-like neutralino is dominated by the fermion channels, one would generally 
expect a relic abundance higher than shown here (where = m^). This is seen 
in Figs. 16 and 17 where an increase in results in larger values of For 

a B-ino-like neutralino, ~ 1 occurs for a neutralino mass of around 100-300 
GeV. 

Along the lower line in Fig. 15 the neutralino is mostly Higgsino. In this 
region, the gauge-boson annihilation channels dominate the total annihilation 
cross section and the results do not depend upon the assumed value of 
cf. Figs. 15-17. However, if the top quark is very heavy (mt w 180 GeV) the 
fermion annihilation channels become comparable to the gauge-boson channels, 
even for a Higgsino-like neutralino. In this case the importance of the fermion 
annihilation channels will depend upon if this ratio is sufficiently large, 

the gauge-boson channels will still dominate. The fact that the Higgsino curves 
in Figs. 15(a), 16, and 17 are the same illustrates this point. In any case, the 
“Higgsino line” in Fig. 17 provides an upper bound to for Higgsino-like 

neutralinos, independent of mj and m^. Note also that our results indicate that 
a Higgsino-like heavy neutralino also provides a viable dark matter candidate for 
a wide range of neutralino masses. 

If < 0.025, it is unlikely that the neutralino is the primary component 
of the dark matter. However, such neutralinos would still reside in our galactic 
halo as a minor component of the dark matter and might still be detectable. 
We should remind the reader that there are few cosmological relics, and even 
relics that contribute only a small fraction of critical density are very interesting. 
Moreover, independent of whether or not the relic neutralinos are the primary 
component of the dark matter, if detected, they still could provide the &st 
evidence for supersymmetry! 
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Unfortunately, a heavier neutralino is more difficult to detect because for 
fixed mass density, the number density is lower for heavier particles. In addi- 
tion, the matrix element for the interaction of a neutralino — or any dark matter 
particle — is typically related to the annihilation matrix element by crossing sym- 
metry, and the annihilation cross section for particle dark matter is always of 
order 10“^^ cm^ sec“^ Since, roughly speaking, the rate for direct detection is 
proportional to the interaction cross section with matter times the relic number 
density and the number density is proportional to the rate for direct de- 

tection is generally inversely proportional to . ’ Likewise, indirect detection 

schemes — searching for the annihilation products of neutralinos that annihilate 
in the halo, in the Sun, or in the Earth — depend upon the annihilation cross 
section times the square of the number density. This implies that for fixed 
the rate for annihilation in the galactic halo falls off roughly as On the 

other hand, the energy of the decay products increases with increasing mj and 
the backgroimds at higher energies tend to be smaller. 

To be more quantitative about indirect detection, in Fig. 18 we show a scatter 
plot of the voliune annihilation rate Fhaio for neutralinos in the galactic halo as a 
fimction of the neutralino mass for a wide range of supersymmetric models. To 
simplify matters, we have assumed that mt = 60 GeV, tan yd = 2, and ttijjo = 0, 
and have used both = rrij^ and mj = oo. Further, we have selected values of 
and M from a grid in the M-fi plane and eliminated models where > 1. For 
the local number density of neutralinos nj, we have used the following simple 
model. For > flhaloi it is reasonable to assume that neutralinos comprise 
the halo dark matter and thvis rij^ — Phalo/”*x- ^ it is very 

unlikely that neutralinos comprise the bulk of the halo dark matter; as a simple 
model we suppose that second case corresponds 

to the situation where neutralinos are a “minor” — ^but potentially interesting — 
component of the halo. To smoothly interpolate between these two cases, we 
have used the simple ansatz. 


PhaJo/^y 


( 7 ) 
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For /Jhalo we use 0.3 GeV cm"^. Based upon the rotation curves of spiral galaxies 
fihalo — 0.03-0.1 (or greater since there is no convincing evidence for the con- 
vergence of the halo mass of any spiral galaxy); for definiteness we have taken 
J^haio^^ = 0 03- The voliune rate for neutralino annihilations in the galactic halo 
is then Fhalo = galactic annihilation rate varies 

roughly as (as argued above). By using m^o = 0, we eliminate the well- 

known s-wave suppression for light neutralino annihilation at zero velocity; 
therefore, the galactic annihilation rates for light neutralinos for larger values 
of mjjo will generally be smaller than those shown in Fig. 18. An interesting 
consequence of our model for galactic neutralino abundance is that the galactic 
annihilation rate does not necessarily decrease with decreasing relic abimdance; 
this is because oc (ctd)”^ and for Z ^^halo, we have fixed rij^. Thus, even 

a neutralino that fails to solve the dark matter problem may prove to be an inter- 
esting candidate for indirect detection. Of course, here we have only considered 
the rate for neutralino annihilation in the galactic halo. To find expected exper- 
imental signatures, the branching ratios for specific final states (e.g., positrons, 
neutrinos, etc.) and the propagation of these annihilation products must also be 
considered.^ 


V. CONCLUSIONS 

In this paper, we have studied in some detail the possibility of a “heavy” neu- 
tralino (m^ > mw), mapping out completely the cosmologically allowed regions 
of parameter space. The main difference between a heavy and a light (mjj < mw) 

neutralino is that additional annihilation channels open up (various pairs of Higgs 

4 5 

and gauge bosons) which have not been included in the previous treatments 
(apart &om a rfecent report of Srednicki and Olive ). We have calculated the 
aimihilation cross sections into these new channels for a general neutralino state. 
In large regions of the parameter spaee for the minimal supersymmetric extension 
of the standard model (/i, M, tan/3, m^o, m^, and mt) the new channels malce 
a substantial contribution to the annihilation cross section, and for a neutralino 
that is Higgsino-like, the gauge-boson final states {W'^W~ , ZZ) often dominate 
all the other diannels by a factor of ten. On the other hand, for a heavy neu- 
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tralino that is gaugino-like, the new final states are typically sub-dominant. The 
new channels contribute and czinnot be ignored for neutralinos that are mixed 
states. 

Using our results for the new annihilation channels along with previous re- 
sults for the fermion final states, we have calculated the total annihilation cross 
section for heavy neutralinos and from this, the relic cosmological abundance of 
heavy neutralinos. We find that for a large portion of the parameter space, a 
heavy neutralino of mass between m\\r and 3200 GeV can have a cosmologically 
interesting relic abundance, t.e., 0.025 < < 1. Thus, we conclude that a 

heavy neutralino is a well motivated and viable dark matter candidate. Based 
upon the cosmological constraint < 1, we have mapped out the regions of 

parameter space that are cosmologically forbidden. In particular, for m< < 180 
GeV, one can completely rule out a neutralino that is heavier than about 3200 
GeV. For a neutralino that is Higgsino-like, the boimd is about 3000 GeV, while 
if the neutralino is gaugino-like, the bound is about 550 GeV. (If a top quark of 
mass less than 120 GeV is discovered, the general boimd drops to 2600 GeV.) It 
is interesting to note that the cosmological upper bound to the neutralino mass 
is comparable to that which follows by insisting that low-energy supersymmetry 
“solve” the hierarchy problem. 
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APPENDIX A: CROSS SECTION FOR xx W+W~ 

First let us review the masses and mixing parameters of the minimal super- 
symmetric standard model. The quantities M, M', and fi are the masses that 
appear in the neutralino mass matrix, Y, Eq. (C38) of Ref. 1, and are the 
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elements of the real orthogonal matrix that diagonalizes Y. The masses of the 
four neutralinos are the absolute values of the eigenvalues of Y, m^o; one should 
keep in mind the fact that the m^o may be negative. The mass of the lightest 
neutralino is the eigenvalue of F which we shall denote by rrij^. The chargino 
masses m-+, given by Eq. (Cl8) in Ref. 1, are always positive. The quemtities <^_ 
and determined from formulas given in Ref. 21, describe the mixing of the 
charginos. The squares of the Higgs-boson masses (m^o and m^o), and 0 are 
given in Appendix A of Ref. 12. Finally, the quantities g and g' are respectively 
the U{1) eind 517(2) gauge couplings. In this Appendix only, some quiintities are 
scaled by mw in order to simplify some of the equations. Specifically, 


w = 




s = 


(j>i +P2? 


m 


W 


(Al) 


where sm^. is the square of the center of mass energy, and p\ and p2 are the four- 
momenta of the incoming neutralinos. (In Appendix B, these same quantities will 
be rescaled by m^.) 

Next we define the quantities 



where e = [ Y| ^ matrix defined in Eq. (C9) of Ref. 1. From these 

we define 

Cf > = 4fl Cf ) = e\U + tifl 

cf) = et + /^ = Cf = e] + fh 

p(^) = ejCj + /i^/l , = eleifi + (A3) 

£>(3) = e^/l -h elfl, = e\e2h + 

= 2eie2/i/2, 

where i = 1,2. 
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The following dimensionless quantities are also useful, 


7 = y(i- 

Li = — log 
75 


Ki = - 


Au) 4 

s s 

(1 + 7 )^ — 2 — 2u? + 2/C| 

_(1 - 7)5 - 2 - 2 u; + 2k,J 

1 


7S [(1 + 7)5 - 2 - 2u; + 2«i (1 - 7)5 - 2 - 2 uj + 2/CiJ ’ 


(A4) 


A = 


me 


sin a + 2[Q sin(/3 — a) — i? sin a] , 

mw 

C = cos a — 2[Q cos{^ — a) + R cos a] , 
tn\\r 


where 


Q=—igZn2-g'Znl), 

9 


R = 


2m w 


{MZl2 + M'Zl,-2flZn3ZnA), 


(A5) 


F = Zl,~ Zl,. 


The total cross section, averaged over neutralino spin for the process XX 
W^W~ is given by 


aww 




Vs-4u>y 


(A6) 


X (Xzz + XffH + ^Bx*) ■ 


The quantities Xu arise from squaring the matrix element given by the Feyn- 
man diag rams in Fig. 2, summing over final state polarizations and averaging 
over initial neutralino spins. Specifically, Xzz comes from the square of the Z~ 
exchange diagram, Xu jj comes from the square of the Higgs exchange diagrams, 
comes from the sum of the squares of the t- and u-channel exchange of 
the charginos, comes from the interference between the t- and «-channel 

exchange of the charginos, and X^-^ and Xg-^. come from the interference of 
the chargino exchange diagrams with the Z and Higgs exchange diagrams. The 
interference between the Z and Higgs exchange diagram vanishes. The quantities 
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Xii axe given by 




Xhh = 


(s - 4)(5 - 4 o;)(s 2 + 205 + 12) 
~2*~ (s - CP + (fz/mz)2 

[2 + (| - 1)“] {» - 4u.)A'", 


8 sin^ /3 


• y fcf >Gf > + cf ^Gf )) , 

X'^X^ 16 ^ * * m\y J 

x~ ~ = ^ £>(^)g('*) + f - i>(^)g(®) (at) 

x+«x+ 8 [ mw yiT^w ^ ^ 

+ !VZ!^d(5)g(7)1, 

m]\r ^w J 




Xrr~. = 


g^F 




S[{s-0^ + i^z/mzyf^^, 


-g^A’ 


8sin^ 

where F 2 ^ is the Z width. Here 
A cos{/3 — a) 


rf^Cp)Gp)+2:^e./,GS'“'), 

" \mw J 

m ' 


A’ = 


Csm(P-a) 


1 1/2 r 1 

1^(5 - m^o/m^)2 + {Tffo/m\vy - ^\olrn^Y + (r£^/mpv)2 


1/2’ 


(A8) 


where F^o is the width of the and the G are given by 


=4KiKi [-(«,• - w - 2)^5 + 4{Ki - 07 - 1)^ + 2u; + 36] 

- i [25^ + 5(40 - 8o7 - 36k.) + 24o7 - 12(«,- - + 120«J 

+ ALi [5 ^ (4k,- — 2k] — K,o7) 

+ s(8 — 0Ki — 12o; + 14 k^ — 6u;«;i + 4o7^ — 3/ci(Ki — o;)^) 

+ 6«, — 8o7 + 9o7^ — 13k] + 6o7k, + 8Ki(K,- — 07)^ — (k,- — 07)^] , 


Gp ^ =24A:i [2 - (k, + ujf - (Ki - 07)2] ^ 16 - 4s 

+ ,-2-t + 2 ^H'“ ' -“)'-«*• + 2^ + ^] 

— 7(/C| — io)^ + 5«i — Tu? + 2| , 
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=2Kil^Ki [(k, — uj — 2)^5 + K, [«i(k» — 4u; — 2) + 6u;^ + 6u; + 5] 

— 2(u> "I" 2)(2w^ — w "I" 2j "I" (cj — l)^(w -j- 2)^ J 

■*" s-2-2uj + 2 k - 8k, • + 4) + 2 s [5«f - (llu) + 14 )k^ 

-|- (7u>^ 16u^ 17 )k,' “■ uj^ — 10u7^ -|- 3tL) — 8j 

+ 4[2 (k, — — 4(k^ — — Kj(K,' — u>)^ 

+ 8k^ 4- 9ojKi + 9k,- — (w — 4)(oJ + 1)] | 

— [— 2s^ + (12 k, — 16oj — 40)s + 36(Kj — — 24 k,- + 88a? + 48] , 

X £i 


= ^ — 2)^K2 ”h (^2 — ”1“ 2(u^ — ^2)^ 

Kj — «2 ^ 

-|- 5^2 — 8^2 — ”H 2)(3c<^ — ^ 

_l_ ? 1^2-s^ + ^ [(^2 ” 6«:2 + 2a? — 3] 

— 4(k 2 — o?)^ + 8«2 + 2a? + 2| + (1 2) 

— [—2s^ + (6/ci + 6«2 ” 16o? — 40)5 
+ 9(^1 + ^2)^ “ (36o? + 12 )(ki + K2) 

+ 3(ki — ^2)^ + 4(9a?^ + 22o? + 12)] , 


G^^^ =8 — 2s — ' [ — (/C2 “ “h 2 — K2 — o?l “h (l ^ 2) 

K\ — K2 

- ®'‘2 - 2“ + 11] 

— 2(/C2 — 4^2 ^ 8o? — 2^ “1“ (l 4— ► 2), 


23 


= (>^2^2 — KiLi) + — — — [6(/tj + K2) + 8u; — 40] 

«i — «2 0 12 *■ 

rl, ^ i2 3. ,2 2 

- [-(ki + K 2 ) - UJ \ + K 1 + K 2 - -(ki - K 2 ) + -W 

+ |«^(4«2 - K2) + ^[-2«2(«2 

2 + 2uj — K\ — K2~ sV 

+ 8/«2 — 8ujk 2 ~ 9^^2 + 4(u? — 2 ){ u > — 1)] 

— (k2 — + 4k 2(«2 — — 5/C2 + (2u; + 2)«2 + 9u;^ — 8a>| +(!<-> 2), 

— [5(^2 — ^ — 2)^ — 4(^2 — 4" 8/C2 — 28w — 4j 4" (1 2) 

/C j *“ ^2 

+ s-2--^+Ki+K2 «2) + » [-(«2 - + 6(<'2 - ‘^)] 

+ 4(«2 ~ “ 8(^2 + 4" 4^ 4" (1 <-» 2) + 2s — 8, 

g|*^ =X, |k,s^(4 + w - K,) 

+ s [4k^ — 5(k, + 3u;) — (k^ — lj)^ — lOw/cj + 6u;^ + 8] 

+ 2 (k, — u;)® — Su}{Ki — u>) — 6Ki — lOw + 4^ 

+ ^ I — s^ + s^(3/c, + w — 18) + s [6(k, — w)^ — 12Kj + 36u; + 28] 

— 12(/C|‘ — — 12ki “I" 44w + 24^, 

=L,|-«,s^ + s[/c, - 2 4- 3w - («j — w)^j 4-4 - 2 (k,- -j-u;) — 2(«:i — u;)^| 

4- — [2s^ -j- 4(K| — w — l)s -f- 8 (k,' — o> -|- 2)] , 

g| ^ = s(l — 5 / 2 ) 4" Li [t(^» "i" w)s^ — s(2 4" 4" 3ti>) 4- 12u?j . (-^9) 


At zero velocity, the Higgs and Z exchange diagrams vanish, and only the 
chargino exchange diagrams remain. The rather lengthy expression for <tww 
times the relative velocity reduces to 


f n1 f 4 + /? , 4 + /I 1^ 


In the early Uixi verse, « 1/4 at freeze-out and the terms are most certainly 
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of importance. However, in the galactic halo, neutralino-neutralino relative ve- 
locities are of order w 0, and to a good approximation, the annihilation 

cross section times relative velocity should be given by (twwv{v -> 0) (unless it 
vanishes). 


APPENDIX B: CROSS SECTION FOR xx 

The matrix element for the process xX is given by the Feynman 

diagrams shown in Fig. 3, and the total cross section, averaged over the neutralino 
spins, is given by an analogous expression to that for aww- 


Xt 

647rsm| 


(r=^) 

In this Appendix, it is most useful to rescale some quantities that have dimensions 
of mass by the mass of the Z. Specifically, 


a; = 




K, 


_ 


mz J \^Z 

where srn^ is the square of the center of mass energy. 
The quantity Xhh is given by 

9^A' ' 




. _ (Pl +P2f 

S — o 5 


ml 


(B2) 


^^^=( 2sinVc ^) 


(B3) 

^2sin^cos^iy/ '■ ' 
where A' is the same as in the W'^W~ case (Appendix A) except that here 
the Higgs masses in the denominators are divided by m| instead of mf^r. (The 
quantities A, C, Q, and R are still given by the same expressions in Appendix 
A.) Next, the quantities Xy'^o and ^x°x°' 

4 


■ -g^A' 

4 sin P cos 6w 




i=l 


me 


Gf)- 

imz rnz 


m^o 


^^(10) 


- l6 ^ ^mz mz ^mz mz mz 


16 ^ 
».j=i 


(B4) 

where the sum is over the four neutralinos. Here the G^j are given by the ex- 
pressions in Appendix A with ki and K 2 being replaced by Ki and Kj (the Li are 
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changed in the same way). The masses of all four neutralinos enter as they are 
virtual particles exch^lnged in the t- and u-channel diagrams. 

The new quantities Ci and Dij that appear above are defined by 


C| — 2e^, Dij — 2e^Cj, (®^) 

where the c, are now given by 

Ci = T — {ZizZnS — Zi^Zni). (B6) 

cos 

For more discussion of the couplings see Refs. 1 and 11 . 

In the zero-velocity limit, the Higgs exchange diagrams again vanish, and 
only the chargino exchange diagrams are nonzero. In this limit 


<^ZZv{v 


0 ) = 


/(U. - 1 )^/^ 
647r|m;f |m2 


E 



1 — w — /c, 


(B7) 


On<% again, this formula is a good approximation for neutralino annihilation in 
the galactic halo, but is not applicable to relic abtmdance calculations. 


APPENDIX C; CROSS SECTION FOR xx 

The total cross section for the process \X HfHj involves the Feynman 
diagrams of Figs. 4 and 5, and is given by 


(TijV = 


9*\k\ 

87rs|mjf| S ’ 


(Cl) 


where |fc| is the magnitude of the outgoing momentum, s = = 4m|/(l— u^/4) 

is the center of mass energy squared, and S is the symmetry factor (2 for i = j 
and 1 otherwise). T^ote that in this Appendix, we have not rescaled any quantities 
by a particle mass. 

The CP quantum number of the final state ^ = c,Cj, where c, is the CP 
quantum number of and ci = C 2 = — C 3 = 1 , determines the form of 
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For ^ = 1 (see Fig. 4), is given by 

= 2p^ 

where 


(B - 2Uf + l(4Ek^Vf 


(C2) 




^ink^jnk 


l=. *0 - 


2m 

+ cim^o + ^777-7-^2 


3(<o - m|o)j 




t=l (to - 

P — 1^(5 — 4 " 


- 1/2 


and Fjyo is the width of Bf. In this case, the contributions from the s-channel 
exchange of B3 and Z vanish. Furthermore, vanishes at r = 0, so that a{jv 
also vanishes as v — > 0. 


For ^ = -1 (see Fig. 5), Xtot is given by 

= —p^k'^{D -^Uf + 2 

3 


m^SD / 

5 \ 

EC I 

1 0 - I - 2V 

2B V 

»T^|/ J 


n 2 


(C4) 


where 


IK 


^ #n - 


^ink^jnk 


Cimy ,8 

" “i£“ 


'« - L 


c = MznnPzijPm^ ^ ~ 


«2 

' + T 


FX 


<2 




2 '\2 io-»T^|o 


u 


2 cos [(s — m|)2 -f m|r|] 


1/2’ 


(C5) 


and is the Z width. In this case, the contributions from s-channel exchange 
of Bj and B® vanish. 
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The expressions given above contain terms of all orders in however, they 
are only correct to 0(v‘^). They are given this way for brevity. By evaluating 
these formulas at two small values of v, the coefficients of the leading order terms, 
a and 6, may be obtained; here, av = a + bv^ + • • 

The kinematic quantities used above are: 


P= l^xb/2, 

= 1 ^ 1 ^ = " 5 + m]jo) + 

to = — ~ 2|mjf |i:, 

t2 = — 2m|, 6 = m^o — m^o. 


The quantities Aij axe the ZHfHj couplings: 


Ail = -^13 = 


sin(a — /5) 
2 cos 6yir 


A32 = —A23 = 


cos(a — 0) 
2 cos 6\y * 


(C6) 


(C7) 


and all other Aij are zero. The quantities Hijt are the couplings, and 

Eijk is symmetric under the interchange of any two indices. The E{jk axe 


o 

Hiu = f— cos(/3 + a) cos 20, 

2 cos u\y 

E 222 = sin(/9 + a) cos 2of, 

2 cos uw 

E\22 = — — [2 sin(y9 + a) sin 2a — cos(y5 + a) cos 2a] , 
2 cos 6 w 

Eu 2 = — ~~§~ [2cos(,S + a) sin 2a + sin(;9 + a)cos2o] , 
2cos»i^ 

cos(^ + a) cos 2^, 

2 COS u\y 

H232 = , + a) cos 2^, 

2 cos 6 w 

E 333 = Eii 3 = E\23 = H 222 — 0. 


(C8) 


The quantities Mjt,> axe the E^xfXj couplings, and Mkij is symmetric in the 
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indices i and j. The M^ij are 


M\ij — 

M20 = 
M30 = 


-1 

'm^oSij sin a 

2 sin j3 

mvy 

-1 

cos a 

2 sin/? 


-1 

cos/9 

2 sin yd 



+ 2{Qij sin(^ - a) — Rij sina')J , 

- 2{Qij cos(/? - o) + Rij cos a)| , 

- 2{Qij cos 2y3 + Rij cos , 


where 

Qij = ^ \Zi^(^gZj2 — g Zji + (i j)] , 

Rij = \MZi 2 Zj 2 + M'ZiiZji — fi{Zi3Zj4 + ZnZji)] . 

2m\\r 

For further discussion of the couplings see Refs. 1 and 11. 

In the limit u — ^ 0 the cross section becomes 


(Tijv{v — > 0 ) 


g^\k\ 


167r|m;f 


C- 


DS 

2m| 


7Tt ry 


) 


‘"-'"a V 4 )\ 


(C9) 


(CIO) 


(Cll) 


APPENDIX D: CROSS SECTION FOR \X ff 

The cross section for neutralino annihilation into quarks or leptons was found 

previously in Refs. 4 and 5, but we include here a slightly improved version. An 

error in the cross section given in Ref. 5 is corrected and Higgs exchange in the 

s-channel is included for the first time. The error involved the neglect of a part 

of the Z propagator and resulted in less than a 10% difference (ignoring Higgs 

exchange) when the cross section was used for annihilation in the early Universe 

(v fa 1/4). However, as pointed out by Lars Bergstrom, it could be significant 

in the v 0 limit which is of relevance for halo annihilation. Higgs exchange is 

relatively unimportant for the heavy neutralino case discussed in this paper, but 

9 

can metke a difference for light neutralinos, especially near poles. 
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To begin, we define the following 


ttq = mgdg/(2Tnw), bg = TziZn2 — i9n6w{TzL — e,)Zni, 

/, = tan^W^e,Z„i, F = Z'^-^ - 

CL = ?3I - Cg sin^ %, CR = -e, sin^ 6wy 

u', = a] + i\, t>J=<^5 + /,^ m’,=a,(h,- f,). 


(Dl) 


Here dg = ZnzJ cos /3 for down- type quarks or charged leptons and dg = Z„ 4 / sin /5 
for up-type quarks or neutral leptons, trig is the mass of the fermion, Tzl is 
the weak isospin, and e, is the fermion electric charge. Leaving out the Higgs 
exchange, the annihilation cross section, averaged over neutralino spins and to 
order is given by 


Cfffv'= ^^^G\cgm\p' [y^ {u'g + v’g){z^-\-{a\ + ri)u^) -j- 4«?'^(l + (02 + r2)t;^) 

9 V '■ 

+ 4w'g{u'^ -f v'^)z(l -I- (-x^ -f- r 3 )u^) + 2u'gVgZ^(l + («4 + 


-h 


7(4 + c%){z^ + ait- 2 ) + ^CLCRz\-l - 04 ^ 2 ) 


+ Fx'^y'^ 


(VgCR - u'gCL){z^ + (ai + z'^rs)v^) 

+ (u'gCR - VgCL)z^{l + (04 + rs)v^) - 2Wg(cL - cr)z{1 + + r5)u^) 




( 




\&m\ 


- 8 -I- 2u' 


2m| 


(1 -I- x"^) - x^ 


-I- F{cl - CR)x^y'^z'^ 2w'gV^ - |ujJt;VV 


ml 


+ 


^2 -1- —2r+ (4u7^ + 2z{u'g + v^)) 


(D2) 
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The quantities ai and rj are given by 


2 b o 1 o 2 ^ /e\ 2 I 2\ 

3 “ 12^ 4^ 02 = -(2-z -i-x ) 

04 = 4" 05 = 1 + 

= I(_4 + ^2 ^ 4 ^ _ ^ ^(_5 ^ 2^2 + 3r^'2 2rz^/3'^) 

3 3 

T3 = ^(-3 + 5r/3'^ - 2^'^) T4 = ^(-3 + 5r^'^) 

rs = ^(-| + r^'^) re = -(-3 - + 2r;5'^), 

(D3) 

where G/- is the Fermi constant, Cg is a color factor (cg = 3 for quarks, Cg = 1 
for leptons), x'^ = m\l{{m\ - sf + r|m|)^/2 is the Z° pole factor, Tz is the 
Z® width, z = mgfmj^, /?' = (! — and y'^ = m^/(m^ + is the 

squark mass suppression including the propagator momentum. The propagator 
momenta factor r = rn^/{m?^ + is usually small, as is — z^). 

When m-j^ — * rrig however, z — > 1, x — ♦ 00 and the expansion breaks down. 
However, this occurs only very close to mass thresholds and so these should be 
avoided. See Ref. 5 for the Feynman diagrams for XX ff further discussion 
of the cross section. 

Besides the squark and Z exchange diagrams included in Eq. (D2), there 
as s-channel diagrams involving the exchange of any of the three neutral Higgs 
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bosons. These contributions can be included by adding the following to Eq. (D2): 
(rfjv = + F2f^'^v^ + 4F3^(1 + 

9 ^ 

2 / 4m| 1,2 4m| \ 

- 2FFzz{cr - cl)x'^ 1 - -i/ + -ri^ - 

+ 2w’,y'^v^ ((F, + F2)f{l - y ) - F,(o 5 + 2re)) 

+ 2z(«; + + Fz)r0\'‘ - 2F3(1 + .-“(rj + j))j 

-Sw'^y'^Fil, 

(D4) 

where Fi = - -s)^ + F^o is the width of Hf, and 


(Qi«,<?2«.,Q3t.) = -^5 (sin a, cos a, cos 

Zm\^ sin /? 

T?T 

(Qu,Q2d,Qid) = ^(cosQ,sinQ,sin/9). 

zm\y cos p 


(D5) 


The quantities Mknn-, tan^, and sin a were defined in Appendixes A and C. 

In the limit u — ♦ 0, relevant for present day neutralino annihilation, the 
annihilation toted cross section (including Higgs exchange) reduces to 


(affv)(v -4 0) = -G%c,ml0' 


- 2 F 3 +y'^ (2u;' + z(u' + u')) 


x'^ 

+ 


(-S) 


(D6) 

For further discussion of the couplings and the supersymmetric model, see Ap>- 
pendix A and Refs. 1, 5 and 11. 
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FIGURE CAPTIONS 


1. Lightest neutralino composition and mass for tan^ = 2. The broken curves 

are contours of constant neutralino mass and the solid curves are con- 
tours of consteuit gaugino fraction -f Z^2)> (a) /z > 0 and in (b) 

/i < 0. 

2. Feynman diagrams for neutralino annihilation into W'^W~ final states. 

3. Feynman diagrams for neutralino annihilation into ZZ final states. The 
index i runs from 1 to 4. 

4. Feynman diagrams for neutralino annihilation into Hj final states where 
Cj = Cj. There are a total of ten diagrams since the index k nms from 1 to 
4. 

5. Feynman diagrzims for neutralino annihilation into HfHj final states where 

a / Cj and only). The index k nms from 1 to 4. 

6. Hatch plots that illustrate the importance of the various channels for neu- 
tralino annihilation in the early Universe {v^ « 1/4). In the areas “hatched” 
with positively sloped lines, annihilation into fermion final states dominates 
by at least a factor of 10; in the regions marked by negatively sloped lines, 
annihilation into gauge-boson final states dominates by at least a factor 
of 10; and in the regions hatched with horizontal lines, annihilation into 
Higgs-boson final states dominates by at least a factor of 10. In the regions 
hatched by vertical lines, several final states contribute comparably to the 
annihilation. Here mj = rtij^ is used throughout. The values of m* and 
tan and the sign of /z are as indicated. Hatch plots for other values of 
mt and tan ^ are qualitatively similar. 

7. Same as Fig. 6(a) but with mj = 2rrij^. 

8. Same as Fig. 6(a) but with = oo. 

9. Relic neutralino abimdance as a ftmction of M for several fixed values of 

/z. The sohd curves show assuming mj — mj, and the broken curves 

show assuming = oo. Here we have taken tan/3 = 2 and mt = 60 
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GeV. Results for // < 0 and other values of tan ^ and are qualitatively 
simil£ir. 

10. Relic neutralino abundance as a function of // for several fixed values of 

M. Ageiin the solid curves show assuming and the broken 

curves show assuming = oo for the same three values of M. Again 
we have taken tan 0 = 2 and mt = 60 GeV. Results for n < 0 and other 
values of tan 0 and m-t are qualitatively similar. 

11. Relic neutralino abundance as a function of M for /x = |M tan^ tan 0 = 

2 and m< = 60 GeV. Again the soUd line shows assuming 

and the broken line shows assuming = oo. Again we have taken 

tan 0 = 2 and mt = 60 GeV. Plots for // < 0 and other values of tan 0 and 
mt etre similar. 

12. .Cosmologically excluded regions of the M-n plane for several values of 

mt using tan^^ = 2. The hatched areas are those for which > 1 

independent of m ffo and m y. Plots for p < 0 and other values of tan 0 are 
very similar. 

13. Same as Fig. 12(a) but assuming mj = oo. 

14. Cosmological upper bound to m^^ as a function of tan 0 for several values 
of the top quark mass. 

15. Scatter plots of vs. neutralino mass for a wide range of models (as- 

smning mj = and tan^ = 2). In (a) we use mt = 60 GeV, and in 
(b) we use m* = 180 GeV. Each “x” represents a different supersymmetric 
model specified by the value of p and M. The values selected for M and fi 
were taken from a grid on the M-fi plane. Results for other values of taii0 
are qualitatively similar. 

16. Same as Fig. 15(a) with mj = 2mj^. 

17. Same as Fig. 15(a) with m^ = oo. 

18. Scatter plot of estimated halo annihilation rate vs. neutralino mass for a 

wide range of models (assuming mjj^ =0, tanj9 = 2, = 60, and using 

= mj or mj = oo). 
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